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We re-examine the question of the entropy stored in the distribution of primordial 
density fluctuations. To this end we make use of two-mode coherent states. These states 
incorporate the isotropy of the distribution as well as the temporal coherence and the 
semi-classical character of highly amplified modes. They also provide a lower bound for 
the entropy if, as one expects, decoherence processes erase the squeezing which originally 
characterized the distribution in inflationary models. This lower bound is one half the 
maximal (thermal) value. By considering backreaction effects, we also provide an upper 
bound for this entropy at the onset of the adiabatic era. 

o 

Q H . I. INTRODUCTION 

o 

Inflation tells us that the primordial density fluctuations arise from the amplification 
of vacuum fluctuations [1, 2]. As a result of this amplification, the initial vacuum state 
becomes a product of highly squeezed two- mode states [3] . In spite of the complexity of 
this state, when computing expectation values, i.e. the two-point function, the modes 
exhibit a temporal coherence upon horizon re-entry. When considering the physics which 
took place near the recombination, it is therefore convenient and sufficient to enforce the 
temporal coherence by putting to zero the decaying mode. Then the residual random 
properties consist in treating the amplitude of the growing mode as a stochastic variable, 
thereby ignoring the quantum properties of the original Gaussian distribution. 

However this simplified description has several drawbacks. In particular the settings 
are too restrictive to describe the distribution which would result from some decoherence 
process which would have taken place during the early universe. More generally, the 
simplified settings are unable to parameterize deviations from the standard results which 
preserve the isotropy and the Gaussianity of the distribution. 

In this paper, we show that the appropriate basis to investigate these questions is 
provided by two- mode coherent states. The reasons are the following. First, a two- mode 
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coherent state provides the quantum counterpart of a particular classical realization of 
the ensemble of metric fluctuations. This correspondence is well defined for the highly- 
excited modes we are dealing with. (Remember that the observed temperature anisotropies 
of relative amplitude 1CP 5 fix the occupation number n to be of the order of 10 100 .) 
Second, ensembles of two-mode coherent states can characterize any Gaussian isotropic 
distribution, thereby allowing to describe arbitrary levels of coherence. This follows from 
the requirement of isotropy and homogeneity which restricts the non-vanishing matrix 
elements of the distribution: only two-mode elements with opposite wave vectors could be 
different from zero. Third, they allow to make contact with the general remark of [4, 5] 
according to which squeezed states rapidly decay into a statistical mixture of coherent 
states when small non-linearities are no longer neglected [6]. In fact when applying this 
general theorem to highly squeezed two-mode states, the resulting distribution is precisely 
a diagonal density matrix of two-mode coherent states. 

This distribution defines the minimal entropy stored in the primordial spectrum (or 
in other words, the minimal coarse graining) when no squeezing remains, that is, when 
there is no direction in phase space in which the spread of the non-diagonal elements of 
the density matrix is smaller than that of the vacuum. For each mode the entropy bound 
is 1/2 the maximal (thermal) value, given the occupation number. We shall see that this 
entropy coincides with that associated with the simplified prescription which consists in 
neglecting the decaying mode. It is to be stressed that this identification can be reached 
because in the large occupation number limit, the entropy is extremely sensitive to the 
level of decoherence whereas the power spectrum is instead extremely robust. Typically 
the relative modifications of the latter are 1/n whereas the changes of the former are in 
Inn. 

The open question concerns the efficiency of decoherence processes in nature: are they 
powerful enough to suppress the squeezing that the initial density matrix possessed? This 
question is currently under investigation. Preliminary results indicate that the squeezing 
is indeed erased, thereby implying that the resulting entropy is larger than (or equal to) 
the above mentionned bound. 

Finally, we shall also provide an upper bound for this entropy by considering backre- 
action effects at the end of the inflationary period. This upper bound is given by 3/4 of 
the maximal entropy. The evaluation of these entropies is exactly performed by exploiting 
the fact that any Gaussian isotropic distribution can be expressed in terms of thermal 
distributions, see Appendix C. 

Related questions have been already analyzed in several papers, see [7-13]. What we 
add in the present paper is a further clarification of the matters, the usefulness and the 
relevance of two-mode coherent states, the lower and upper bounds on the entropy, and 
relationships between various elements which have been some how separately discussed. 
Notice that we shall not discuss the physical relevance of this entropy for structure forma- 
tion. For this interesting question we refer to [7]. 
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II. A REVIEW OF THE STANDARD DERIVATION OF PRIMORDIAL 

SPECTRA 

A. Quantum distribution of two-mode states 

In this subsection we recall how the amplification of modes of a massless field propa- 
gating in a FRW universe translates in quantum settings in the fact that the initial ground 
state evolves into a product of highly squeezed two-mode states. Before proceeding, we 
remind the reader that it has been shown that the evolution of linearized cosmological 
perturbations (metric and density perturbations) reduces to the propagation of massless 
scalar fields in a FRW spacetime [14]. In this article, we shall only consider the mass- 
less scalar test field since the transposition of the results to physical fields represents no 
difficulty. Indeed, when preserving the linearity of the evolution, the only modification 
concerns the late time dependence of the modes. 

Let us work in a fiat FRW universe. The line element is: 

ds 2 = a{r]) 2 [-drj 2 + SijCbfdat] . (1) 

For definiteness and simplicity, we consider a cosmological evolution which starts with an 
inflationary de Sitter phase and ends with a radiation dominated period. When using the 
conformal time r\ to parametrize the evolution, the scale factor is respectively given by 

a{r)) = - H(r,-2r, r ) > ^ " °° < V < Vr ' (2&) 
a ^ = JT^ V ' for T l>Vr, (2b) 

±1 Tj r 

where r\ r > designates the end of inflation. The transition is such that the scale factor 
and the Hubble parameter are continuous functions. This approximation based on an 
instantaneous transition is perfectly justified for modes relevant for CMB physics. Indeed, 
their wave vector k obeys kr\ r ~ 10~ 25 ~ e" 60 when inflation lasts for 60 e-folds. Hence, the 
phase shift they could accumulate in a more realistic smoothed out transition is completely 
negligible. 

Let £(77, x) be a massless test scalar field propagating in this background metric. It is 
convenient to work with the rescaled field 4> = a£ and to decompose it into Fourier modes 

<Kr?,x) = Jd 3 k^^^( V ). (3) 



The time dependent mode 0k obeys 

^+(^-^14 = 0, (4) 



dla 



where k = |k| is the norm of the conformal wave vector. 

In our background solution, k 2 — d 2 a/a is negative during the de Sitter period when 
the wavelength is larger that the Hubble radius. This leads to a large amplification of 
In quantum terms this mode amplification translates into spontaneous pair production 
characterized by correspondingly large occupation numbers. 
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To obtain the final distribution of particles, one should introduce two sets of positive 
frequency solutions of Eq. (4). The in modes are defined at asymptotic early time, and 
the out ones at late time. Both have unit positive Wronskian in conformity with the usual 
particle interpretation [15]. One gets 

whv) ■ for " < * • <5a> 

C'(>ri = 4r^ i "'. for (») 

\J2k 



In spite of the time dependence of the background, these two positive frequency modes are 
unambiguously defined (up to an arbitrary constant phase which drops in all expectation 
values and which has here been chosen so as to simplify the forthcoming expressions). In 
the radiation dominated era there is no ambiguity since the conformal frequency is constant 
because d^a = 0. In the de Sitter epoch, there is no ambiguity either for relevant modes if 
inflation lasts more than 70 e-folds, see [16] for the evaluation of the small corrections one 
obtains when imposing vacuum at some finite early time. Similarly, in quantum settings, 
there is no ambiguity for the initial state of relevant modes: at the onset of inflation they 
must be in their ground state[14, 17]. 

The in and out modes are related by a Bogoliubov transformation 

W{ V ) = a k 4F{v) + ft<ft* m {v), (6) 
where the Bogoliubov coefficients are given by the Wronskians 

These overlaps should be evaluated at transition time r\ r since modes satisfy different 
equations in each era. One gets 

2ikrj r _ ^ 

= -2khi (1 " 2ikn ' ~ = 2kh§ (1 + ° {kn ' )3) • <8a> 

(8b) 

Thus, for relevant modes, krj r ~ 10 -25 , the in modes are enormously amplified. Concomi- 
tantly, they are dominated by the sine during the radiation dominated era: 



sin kf] 



+ 0(krj r ) cos krj 



V>Vr- (9) 



Once the cosine in Eq. (9) is neglected, the physical modes <^l n /a show the same temporal 
behaviour, e.g. they are constant until they start oscillating as they re-enter the Hubble 
radius when kr) ~ 1. 

Lets now see how these considerations translate in second quantized settings. Each 
mode operator is decomposed twice 

Mv) = K4(v) + ^A*(v), (io) 
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where j stands for both the in and out basis. The operators so defined are related by the 
transformation 

a™ = a* k a^-(3 k a^. (11) 

This transformation couples k to — k only. Hence, when starting from the in vacuum (the 
state annihilated by the operators), every out particle of momentum k will be accompa- 
nied by a partner of momentum — k. Moreover, pairs characterized by different momenta 
are incoherent (in the sense that in expectation values any product of annihilation and 
creation operators of different momenta will factorize). 

These two properties are explicit when writing the in vacuum in terms of out states 
(i.e. states with a definite out particle content). From Eq. (11), one gets (see [7, 19]) 

|0, in) = Y[ ® |0, k,m) 2 

k 

= n ( ex p { zk a k u<t &ou k) i°> k > ° ut ) ® i°' - k > ° ut )) ■ ( i2 ) 

k 

The tilde tensorial product takes into account only half the modes. It must be introduced 
because the squeezing operator acts both on the k and the — k sectors. The definition 
of this product requires the introduction of an arbitrary wave vector to divide_the modes 
into two sets. The sign of k z can be used. Notice that a rigorous definition of J^k requires 
to consider a discrete set of modes normalized with Kronekers (that is, to normalize the 
modes in a finite conformal 3-volume). To be explicit, the two-mode state |0, k, in)2 is 
given by 

|0, k, in) 2 = |0, k, in) ® |0, -k, in) , (13) 

where |0, k, in) is the ground state of the k-th mode at the onset of inflation. The complex 
parameter z k appearing in the squeezing operator in Eq. (12) is given by the ratio of the 
Bogoliubov coefficients 

z . - Ph. = _ e -2ifcr? r 1 

a* " (l + i2k7] r -2k 2 rft) 

= -l + 0(k Vr ) 3 . (14) 

The high occupation number limit corresponds to \zk\ — > 1~. 

It has to be emphasized that none of the out states in Eq. (12) carry 3-momentum. 
Hence, the distribution is homogeneous in a strong sense: at late time the 3-momentum 
operator is still annihilated by the state of Eq. (12). (This property is not satisfied by 
incoherent distributions such as thermal baths. In those cases, the 3-momentum fluctuates 
and vanishes only in the mean.) The present distribution is also isotropic since the Bogoli- 
ubov coefficients are functions of the norm k only. Finally, it is a Gaussian distribution, 
as can be seen from Eq. (12). 

To appreciate the peculiar properties of the distribution of Eq. (12) it is interesting 
to consider the most general homogeneous, isotropic, and Gaussian distribution of out 
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quanta. Its properties are completely specified by three real functions of the norm k (one 
real and one complex) through the following expectation values 

<ar*>=0, (15a) 
(a™ tt aff t ) = n k S 3 (k-k'), (15b) 
= c fe <5 3 (k + k'). (15c) 

In the second line, n k is the mean occupation number. In the third one, the complex 
number c& characterizes the quantum coherence of the distribution. The degree of two- 
mode coherence is given by \c k \/{n k + 1/2), see Appendix C. It is bounded by 1. For a 
thermal (incoherent) distribution, one has Ck = 0: no coherence. 
In the case of pair production from vacuum, one has 

I |2 
I a 1 2 \ z k\ 

nk — \Pk\ —~ : — [2 , 

1 - \Zk\ 

c k = a k (3 k = • (16) 

1 - \Zk\ 

Therefore, when considering relevant modes in our inflationary model Eq. (8), one has 

„. = I ~ io 100 

f;l - l + 0(k Vr f. (17) 



n k + 1/2 



That is, the distribution which results from inflation is highly populated and, more impor- 
tantly, maximally coherent. When computing the Green function, the two-mode coherence 
of the distribution will manifest itself in the temporal coherence of the modes. 



B. Two-point function and the neglect of the decaying mode 

When expressed in terms of out modes, the two-point function associated with the 
general distribution specified by Eqs. (15) is 

G(n,x, r/,x' ) = G out (n, x, n' , x') 

+ r Sin , ( f |X " ^ 2Re [n k Mv)<f>Uv') + c k Mv)Mv')] ■ (18) 

JO TT K\X — X | 

In the first line, we have isolated G ou t, the contribution of the out vacuum. In the high 
occupation number limit, this contribution is negligible. 

It is important to notice that, for a general distribution, the sum in bracket in the 
above integrand cannot be factorized. However it does factorize in two cases which are 
relevant for us: first, for coherent states, see Appendix A, and second for distributions 
resulting from pair production from vacuum. Indeed, when considering our cosmological 
model, taking into account the minus sign of Eq. (14), and neglecting correction terms in 
0{krj r ) 3 (which amounts to neglect the decaying mode, see Eq. (9)), one obtains 

_ , , f°° dkk 2 sin(fc|x — x'l) sinfcr? s'mkri' 
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As announced, the integrand in Eq. (18) factorizes. In inflation, the function which 
appears is sin krj where 77 is related to the scale factor by Eq. (2). This is how the 
temporal coherence of modes obtains from the two-mode coherence of the distribution. 

Once the cosine is neglected, the quantum distribution is effectively replaced by a 
stochastic Gaussian distribution of classical fluctuations 

, . . „ sin/cr? , , 

Mv) = s ^-^ > ( 20 ) 

with locked temporal argument, and random amplitudes with variances given by 

(( S k S£, » e// = (( S k S„ w )) eff = 2n k 6® (k - k') . (21) 
Being Gaussian, the effective probability distribution is simply 

*«=fi£«p(-tr)- < 22 > 

To avoid double counting, one must again use the tilde product which takes into account 
half the modes only, as was done in the quantum distribution of Eq. (12). This counting 
becomes crucial when computing the entropy, see Section III.C. 

Notice finally that the first equality in Eq. (21) is simply the expression of the reality 
of the field (p(rj, x). This was not the case when noticing that inflation gives (d£ u ^ a™*) = 
— (a£ u * a°^/)(l + 0(kr] r ) s ) for the quantum field operators. In that case the equality is the 
expression of the coherence of the in vacuum, i.e. the absence of 3-momentum fluctuations 
and the possibility of factorizing the 2-point function by discarding the cosines. 



C. Additional remarks 

First we remind the reader why a random distribution of both the sine and the cosine 
does not give rise to any temporal coherence [3]. In fact such a distribution corresponds 
to an incoherent (thermal) distribution. 

Writing the field modes as the sum of a sine and a cosine: 

~ cos(fcr/) q sin(fcr?) 

= Cl wr + Sk ^/r ■ (23) 

It is to be noticed that the operators Ck and Sk are proportional to the field mode operator 
and its conjugate momentum evaluated at r] = 0, as if there were no inflation: 

C k = (o^ + alf) =VkM0), 

5 k = (ar - air) = 7^ d v M0) • (24) 



They satisfy the canonical equal time commutation relations. 
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Consider an incoherent {c k = 0) distribution: 

(ala^)^ = n k 5 3 (k - k') , (a^a^)^ = . (25) 
Then Ck and are uncorrelated Gaussian operators with equal variance: 

(C k< S k '>^ = 0. (26) 

This implies the absence of temporal coherence of the modes, as can be seen from the 
temporal behaviour of the bracket in the integrand of Eq. (18): when r] = rf the bracket 
does not exhibit any oscillation in k as it did for the inflationary distribution. (This absence 
can also be understood by considering the 5k and Ck as stochastic variables rather than 
quantum ones. Writing the mode in terms of its norm and its phase [3, 20] 

fa = $ k sm(krj + k ) . (27) 

one can treat <l>k and #k as stochastic variables. One then verifies that the distribution 
for the phase k is uniform over the interval [0, 2tt]. Hence no temporal coherence could 
obtain.) 

As a last remark, to estimate what has been neglected when discarding the cosine in 
Eq. (20), we compute the residual fluctuations of Ck and the cross correlation C^S^. To 
leading order in n k , one has 

{SJl) in = (n k + X - - Re(c k )\ 5 3 (k - k') = 2n k 5 3 (k - k') , 

(C k Cl) m = (n k+ l - + Re{c k ) s j 5 3 (k-kO = 0(ni/V(k-k'), 

(C k Sl) m = - (lm{c k ) + ^ 5 3 (k - k') = 0(nl /4 ) 5 3 (k - k') . (28) 

When divided by the variance of 5k, the variance of Ck and the cross-correlation are of 
order n k 3 ^ ~ {kr) r ) 3 ~ 10^ 75 . One can therefore safely use the distribution Eq. (22) in 
replacement of the quantum distribution Eqs. (12) or (28) when calculating the power 
spectrum. However, in general, this is not the case for the entropy. 



D. Drawbacks of the simplified description 

The simplified description in terms of a statistical ensemble of sine standing waves has 
indeed several shortcomings. It is of value to describe them with some attention. 

First, it should be pointed out that in the early universe, different physical processes 
could give rise to different final density matrices, and hence, to different entropies. As 
we shall see bellow, only a narrow set of these quantum distributions can be put in cor- 
respondence with that of Eq. (22). To describe the most general isotropic and Gaussian 
distribution, it is necessary to return to two-mode distributions characterized by n k and c k 
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of Eq. (15). Second, the classical ensemble of sine waves should be considered only as an 
effective description of the density matrix. The appropriate basis to describe this matrix 
is provided by coherent states. The reasons are the following. On the one hand, they 
constitute the quantum counterparts of classical configurations in phase space. Therefore 
they provide an adequate basis for studying the semi-classical limit. In particular, since 
the spread of coherent states have no preferred direction in phase space, it will be easy to 
see whether or not a distribution has kept some squeezing. On the other hand, they are 
the preferred basis in which squeezed states decohere when weak interactions are taken 
into account [5]. 



When using coherent states in cosmology, one must pay attention to the entanglement 
between k and — k modes. A naive use of coherent states which would assign amplitudes to 
each mode separately could erase these correlations and therefore suppress the information 
about the temporal phase of the modes. Taking into account the entanglement leads to 
the notion of 'two-mode coherent states'. 

A. Representation of the in vacuum with two-mode coherent states 

To understand the usefulness of two-mode coherent states it is appropriate to first men- 
tion the following properties [21]. Consider a mode k in a coherent state \v, k) constructed 
with out operators, see Appendix A for its definition. Then compute the one-mode reduced 
state obtained by projecting it on the two- mode initial vacuum of Eq. (13) 



It is remarkable that the state of the — k mode which is entangled to \v, k) is also a coherent 
state. Its amplitude is given by the complex conjugate of v times Zk characterizing the 
pair creation process. These facts follow from the EPR correlations in the initial vacuum 
displayed in Eq. (12). The prefactor Ak(v) is 



It is the probability amplitude to find the mode k in the coherent state \v , k) given that 
we start with vacuum at the onset of inflation. 

Using the representation of the identity with coherent states, Eq. (A13), the two-mode 
in vacuum can be thus decomposed as a single sum of two-mode coherent states rather 
than two independent integrations over one-mode coherent states. In fact we have 



III. TWO-MODE COHERENT STATES 



{v, k|0m, k) 2 = Ak(v) \z k v, -k) . 



(29) 




(30) 




(31) 
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This result is exact and applies even for low occupation numbers. It is the consequence 
of the coherence of the in vacuum and holds for every homogeneous pair creation process. 

Another important consequence of Eq. (29) is that the probability to find simultane- 
ously the k-mode with coherent amplitude v and its partner with amplitude w is 

V 2 ,k(v,w) = \(v, k\(w, -k|0m, k) 2 | 2 = \A k (v)\ 2 x e~\ w ~ z ^ . (32) 

The second factor arises follows from the overlap between two different coherent states: 
|(u|f)| 2 = exp(— \u — v\ 2 ). Equation (32) implies that once the amplitude of the k-mode 
has been measured, the conditional probability to find its partner in a coherent state \w) 
is centered around w = z k v* . In the high occupation number limit we are dealing with, 
the spread (= 1) around this mean value is negligible when compared to the spread in v 
which is given by \ak\ = n& + 1. Therefore, when computing expectation values in leading 
order in n k , the conditional probability acts as a delta function on both the real and the 
imaginary part of w. This is how the EPR correlations in the in- vacuum translate in the 
coherent states basis. This result determines the properties of the local correlations in the 
primordial spectra[21]. 

To complete this analysis, and in preparation for studying decoherence, it is also in- 
teresting to explicitly write the non-diagonal matrix elements of the in vacuum density 
matrix. One has 

(v\(w\p in \v'}\w'} = A2,k(v,w)A 2 ,k{v',w')* , (33) 
where the two-mode amplitude is 

A 2 ,k(v,w) = A k (v)e-^ w - z * v *\ 2 e ilm(w * z * v *K (34) 

Since the initial vacuum is a pure state and since [a^l 3> 1, the above matrix elements do 
not vanish, even for macroscopically different coherent states. Therefore this distribution 
does not describe a classical ensemble of these quasi-classical states. Fortunately, such 
quantum distributions are unstable to any weak perturbation in that they rapidly evolve 
into statistical mixtures. Let us now describe this decoherence process. 

B. Zurek et al. analysis and minimal decoherence scheme 

In general, it is a difficult question to determine into what mixture an initial density 
matrix will evolve when taking into account some interactions amongst modes or with 
other modes. There exist however several cases where clear conclusions can be drawn. 
First, when one can neglect the free Hamiltonian, the preferred states (that is the basis 
into which the reduced density matrix will become diagonal) are the eigenstates of the 
interaction Hamiltonian [22-24]. This approach has been applied in [12], to primordial 
density fluctuations when the (physical) modes are almost constant because their wave 



1 Notice however that the above decomposition is not unique since the coherent states are not orthogonal, 
compare with [27]. Eq. (31) has the advantage to be directly related to the detection of a quasi-classical 
configuration in the k sector. 
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length is much larger than the Hubble radius. The conclusion is that the preferred ba- 
sis is provided by amplitude (position) eigenstates. However this conclusion leaves some 
ambiguity and might lead to some difficulties. First, position eigenstates are not nor- 
malized. Second, and more importantly, the spread in momentum is infinite for these 
states. Therefore, the velocity field would not be well defined when the modes re-enter the 
horizon. Moreover, as pointed out in [13], some additional decoherence could be obtained 
as they re-enter the horizon. In this case, the momentum should be treated in the same 
footing as the position. To cure these problems, some finite spread in position should be 
introduced. One then needs a physical criterion to choose this spread. 

To remove this ambiguity it is appropriate to appeal to coherent states both for mathe- 
matical and physical reasons. In this article, we shall only present the basic mathematical 
results. We reserve for a forthcoming publication the justification of the physical rele- 
vance of these states in a cosmological context. Let us simply notice the following points. 
In inflationary cosmology, modes are weakly interacting harmonic oscillators [6]. Indeed, 
given that the relative density fluctuations have small amplitude (~ 10~ 5 ), the hypothe- 
sis of weak interactions is perfectly legitimate. Second, coherent states provide the basis 
in which the density matrix decoheres when considering weakly interacting harmonic os- 
cillators. This has been shown by Zurek & al. [5]. The criterion they used to reach 
this conclusion is the minimization of the growth of entropy in the course of the evolution. 
With this criterion, coherent states are more stable than squeezed states in that the growth 
of entropy one obtains when they are used as initial states is much slower. Hence, when 
starting with a squeezed state, there is a phase of rapid growth of the entropy which sends 
the system into a mixture of coherent states and which is followed by a period of slower 
increase. The entropy growth is in fact directly related to the decay of the squeezing. 
Since it is unlikely that the interactions in the early cosmology could be sufficiently weak 
so as to keep some squeezing, we can use the following mathematical result to infer that 
the actual entropy of the final distribution should be higher than (or equal to) a certain 
bound. 

Coherent states indeed define a minimal decoherence scheme in the following sense. 
Consider the set of final distributions which result from the initial distribution of Eq. (33) 
through some decoherence process and which no longer possess any squeezed direction. 
The lowest value of the entropy in this set is given by the entropy of the incoherent su- 
perposition of coherent states, with statistical weights given by the probabilities to find 
the corresponding coherent states, as in Eq. (32). This distribution gives a lowest entropy 
simply because coherent states have minimal constant spreads (given quantum uncertain- 
ties and when considering free evolution, see Appendix A). We shall now explicitely write 
down this distribution and compute the entropy it carries. 

C. Application to cosmology 

When considering the initial distribution Eq. (33), one obtains the following density 
matrix: 

/dPv dPw 
V 2 ,k(v,w) \v, k)(v, k| (g> \ w, -k)(u>, -k| , (35) 
TT IT 
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where the probability distribution is given in Eq. (32). In Appendix B we show that 
Eq. (35) is indeed the resulting normalized distribution. The technical point which re- 
quires clarification is the extension of the above mentionned minimal scheme to two-mode 
squeezed states. 

It should first be noted that when computing expectation values in leading order in n k , 
the above distribution can be simplified and written as a single sum of two-mode coherent 
states, as in eq. (31): 

\A k (v)\ 2 \v, k){v, k| ® \zv*, -k){zv*, -k| . (36) 

As we shall progressively see, this distribution should be conceived as the quantum coun- 
terpart of the effective distribution of sine functions discussed in section II. B. 

Secondly by reducing the density matrix, some entropy has been introduced. One 
verifies indeed that Tr(p^ iri ) < 1. The important point is that this decohered distribution 
has retained all the information about the temporal coherence of the modes. In fact one 
has 



Zk\oik\ 2 ■ 



m / - - - \ fdvdw . . , fd v r 2 , * n 

Tr(p min a k a_ k ) = / [P 2 , k{v, w) v w\ = / \A k {v)\ v [z k v , 

J IT IT J TT l 

t t 2 

Ti(p min ala k ) = Tr(p min a]_ k a_ k ) = \a k \ =n k + l. (37) 

The first line shows that the cross term is equal to that of the original distribution, see 
Eqs. (15, 16). In the second line, one sees that the occupation numbers slightly differ, but 
to order l/n k only. 

The effect of having increased by 1 the occupation number while having kept untouched 
the cross-correlation means two things. First, the (relative) degree of coherence has been 
reduced and therefore some entropy has been created. Secondly, in the high occupation 
limit, the two-point function of Eq. (18) is not affected by this loss of coherence: for 
relevant modes, the relative change being of the order of l/n k ~ 10~ 100 . 

D. Minimal entropy and the neglect of the decaying mode 

The entropy of any Gaussian two-mode distribution can be exactly calculated [25, 26] by 
using the fact that the density matrix of a two-mode squeezed state is unitarily equivalent 
to the tensorial product of two thermal density matrices of oscillators, see Appendix C. 
We shall name a and b these two real oscillators. One has 



P2,k 



M* pth,a®Pth,bM, (38) 



where M is a unitary operator acting on the two-mode Hilbert space. The expression of 
the (von Neumann) entropy immediately follows: 

S[h,k] = S[p tKa } + S[p tKb ] , (39) 

where the entropy of a thermal bath with mean occupation n is 

S [pth] = (n + 1) ln(n + 1) - ft ln(n) . (40) 
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When considering only distributions preserving homogeneity and isotropy, the occupation 
numbers of the thermal matrices are equal and given by 

If 1 \ 1/2 

n k + ^ = Un k + -) 2 - \c k \ 2 J , (41) 

where n k and c k are defined in Eq. (15). 

Let us apply this result to several cases. First, for the two-mode in vacuum of Eq. (12), 
the occupation number and the coherence term are related by Eq. (16), one has n k = 
as expected. Hence the entropy vanishes. 

For the decohered matrix Eq. (35), using Eq. (37), the occupation number of the two 
thermal baths are 

n k = g (- 1 + V8(n fe + 1) + l) ~ ^hT k , (42) 

where the last term is the leading order when n k S> 1. The two- mode entropy of this 
mixture is then 

S[pmin] = 2 5 [pth] = 2\nn k + O(l) , 
= \nn k + 0{l) =2r fc + 0(l), 

~ 1001n(10). (43) 

In the second line, we have expressed the occupation number in term of the squeezing 
parameter r k : n k = sh 2 r k . Hence, a two-mode squeezed vacuum state which decoheres 
in the two-mode coherent basis goes along with an entropy of 5k, -k = 2 r k per two-mode. 
This value is large, but not maximal. Indeed, had the coherence term c k vanished while 
keeping the occupation numbers fixed [7], one would have found the maximal value of the 
entropy which is given by twice this above value, i.e. 

^ k = S£i k = 4r fc , (44) 

or 5 k nc = 2 r k per mode k. 

It is interesting to notice that the entropy associated with the Gaussian distribution 
Eq. (22) of sine functions equals that of Eq. (35), up to an arbitrary constant which 
arises from the usual ambiguity of attributing an entropy to a classical distribution. (This 
ambiguity can be lifted when introducing h to normalize the phase space integral.) Using 
this trick, the entropy associated with Eq. (22) is 5 e // = 2 r k for each independent mode, 
since this entropy is maximal. However, for each independent mode here means for each 
two-mode since the mode — k is no longer independent of the k mode once the cosines 
have been neglected, see Eq. (21). From this equality of entropies, we conclude that the 
quantum density matrix which corresponds to the Gaussian ensemble of sine functions is 
precisely given by Eq. (35). 

A priory one might think that many quantum distributions can be associated with the 
classical distribution Eq. (22). This is not the case when imposing that Gaussianity is 
preserved and that the entropies coincide. Indeed, in the high squeezing limit, the entropy 
is an extremely sensitive function of the relative coherence. To see this dependence, let 
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us calculate the entropy of a generic distribution (15), and let us write the norm of the 
coherence term as 

\ck\ 2 = n k {n k + 1 - 5 k ) , (45) 

where 5 k is a real number between and n k + 1. The equation (41) has the solution 
n(n + 1) = n k 5k- They are three characteristic values. 5k = obviously corresponds to 
the pure squeezed state: the in vacuum. 5k = 1 corresponds to the minimal decoherence 
scheme with entropy given in Eq. (43). 5 k = + 1 corresponds to the thermal case 
with maximum entropy. From this analysis we see that the uncertainty in the definition 
of the quantum distributions which give rise to the entropy 5k, -k = 2r^ is very limited: 
5^ must be of order 1 in the following sense. Consider that the loss of coherence scales 
as 5k oc nj.. Then the thermal occupation number and the entropy respectively scale as 

fik oc n^ +1 ^ 2 and Sk, -k = (1 + 7)2r^ + Const. This linear dependence in implies that 
the distributions with entropy given by Eq. (43) all have 7 = 0. 

We re-emphasize that a value of 5k smaller than 1 is unlikely in the context of primor- 
dial fluctuations since it would mean that the distribution has kept some of its quantum 
squeezeness. The remaining question thus concerns the computation of 5k, noticing that it 
can receive contributions both from the inflationary period and from the adiabatic era[13]. 
The challenge is to determine which one is more important and what could be a realistic 
value of 5k- 

As a final comment, we provide an upper bound for the decoherence entropy which could 
have resulted from processes in the inflationary phase. Because increasing the decoherence 
implies increasing the power of the growing mode, one obtains a bound on the decoherence 
level when requiring that the power of the decaying mode be equal to that of the growing 
mode at the onset of the adiabatic era. (This requirement follows from the fact that 
the rms value of the primordial fluctuations (of the Bardeen potential) cannot be much 
higher than that obtained from in vacuum because otherwise this would invalidate the 
whole framework of linear metric perturbations.) Using Eq. (18) evaluated at r\ r and the 
parameterization of Eq. (45), one obtains 

5k = nf , S^_Z = 3r k + O(l) . (46) 

If no further decoherence is added in the adiabatic area, this should be the maximum 
amount of entropy stored in the primordial spectrum. Notice that when evaluated at 
recombination, the two-point function is still unaffected by this modification of the co- 
herence because at that time the decaying mode has still further decreased. Indeed the 
residual modifications are then of the order of n^ 1 ^ 2 ~ 1CT 50 . 

Acknowledgments: We are grateful to Dani Arteaga, Claus Kiefer, Serge Massar, 
Jihad Mourad and Alexei Starobinsky for useful remarks. 

APPENDIX A: COHERENT STATES 

This appendix aims to present the properties which we shall use in the body of the 
manuscript. For more details, we refer to [28-30]. 
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Coherent states (of a real oscillator) can be defined as eigenstates of the annihilation 
operator: 

a\v) = v\v) , (Al) 
where v is a complex number. In Fock basis it is written as 

H = e ~^E4ii n >' ( A2 ) 



n=0 v n 

where the exponential prefactor guarantees that the state is normalized to unity {v \v ) = 1. 
They are also obtained by the action of a displacement operator on the vacuum : 

\v) =D(v)\0) = e "*" t . (A3) 

The first interesting property of coherent states is that they correspond to states with 
a well defined complex amplitude v. Indeed, by definition (Al), the expectation values of 
the annihilation and creation operators are 

(v\a\v) = v , {v\a i \v)=v*. (A4) 

Thus the mean occupation number is 

(v\a^a\v) = \v\ 2 . (A5) 

It is to be also stressed that the variances vanish: 

Aa 2 = (v\a 2 \v) - (v\a\v) 2 = 0, Aa f2 = {v\a) 2 \v) - {v\a*\v) 2 = . (A6) 

From these properties one sees that the expectation values of the position and momen- 
tum operators (in the Heisenberg picture) 

q(t) = (ae~ iujt + aV^) , p(t) = -i^J^(ae~ iu}t - a)e iujt ) , 



are 



q(t) = {v\q(t)\v) = J^(ve~ iujt +v*e iujt ) = J—\v\cos(cot 
V 2uj V ui 




pit) = (v\p{t)\v) = -i\ — {ye~ iwt - v*e iut ) = -y/2hu\v\ sin(wt - <p v ) = d t q(t) .(A7) 



We have used the polar decomposition v = \v\e l< ^ v . These expectation values have a 
well defined amplitude and phase and follow a classical trajectory of the oscillator. This 
is due to the "stability" of coherent states under the evolution of the free Hamiltonian 
2Hq = p 2 + uj 2 q 2 : if the state is \v) at time to, one immediately gets from (A2) that at a 
later time t, the state is given by \v(t)) = ^e^ 1 ^^ 1 ^) . Notice that the variances of the 
position and the momentum are 

Af = -, Af = -. (A8) 
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They minimize the Heisenberg uncertainty relations and are time-independent. Hence, in 
the phase space (q,p), a coherent state can be considered as a unit quantum cell 2-kH in 
physical units (see also (A14) for the measure of integration over phase space) centered on 
the classical position and momentum of the harmonic oscillator (q(t),p(t)). In the large 
occupation number limit \v\ S> 1, coherent states can therefore be interpreted as classical 
states since Aq/^q 2 + p 2 /tu 2 = Ap/^to 2 q 2 + p 2 = l/2\v\. This is a special application 
of the fact that coherent states can in general be used to define the classical limit of a 
quantum theory, see [30] and references therein. 

One advantage of coherent states [28] is that the calculations of Green functions re- 
sembles closely to those of the corresponding classical theory (i.e. treating the fields not 
as operators but as c-numbers) provided either one uses normal ordering, or one considers 
only the dominant contribution when \ v\ 3> 1. We compute the Wightman function in the 
coherent state \v) 

G v (tX) = (v\q(t)q(t')\v) 

= (: mW) :>„ + 1;^-'') , (A9) 

where we have isolated the contribution of the vacuum. The normal ordered correlator is 
order \v\ 2 : 

(■ q{tW) -)v = -Re \{a 2 ) v e~ w ^ + {a)a) v e w ^ 

UJ L 

= — \v\ 2 cos(u)t - <f> v ) cos(ujt' - <f) v ) = q(t)q(t') . (A10) 

We see that the perfect coherence of the state, namely |(aa)„| = {a)a) v is necessary to 
combine the contributions of the diagonal and the interfering term so as to bring the 
time-dependent classical position q(t) in Eq. (A10). 

The wave-function of a coherent state in the coordinate representation is given by 

where v = (ujq + ip) / V 2ujH. This follows from the definition (q\a\v) = v(q\v). From this 
equation one notes that two coherent states are not orthogonal. The overlap between two 
coherent states is 

(v\w) = exp ^v*w — t^M 2 — ^l^l 2 ^) ' (A12) 

Nevertheless they form an (over)complete basis of the Hilbert space in that the identity 
operator in the coherent state representation {\v}} reads 

1 



h-^\v){v\. (A13) 



The measure is 

d 2 v d(Rev)d(lmv) dqdp 



(A14) 



7r 7r 2irh 

The representation of identity can be established by calculating the matrix elements of 
both sides of the equality in the coordinate representation {\q}}, with the help of (All). 
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APPENDIX B: APPLICATION OF ZUREK & AL. RESULTS TO THE 

COSMOLOGICAL CASE 



In this appendix we show that Eq. (33) is indeed the minimal decohered distribution 
by decomposing the complex mode <^>k into two real oscillators <pi and <f>2 given by its 
real and imaginary parts. Since the two-mode Hamiltonian is Hermitian, it splits into the 
sum of two identical one-mode oscillator Hamiltonians for (pi and 02 separately. Notice 
that the annihilation operators of these two real oscillators, a\ = (ctk + fi-k) /V%, «2 = 
—i (at — «-k) /V% rnix k and — k annihilation operators. Hence they can easily take into 
account the entanglement between these two modes. 

A two-mode squeezed state |0m,k)2 can always be written as the tensorial product of 
the two one-mode squeezed states[31]. In our case, the one-mode squeezed states are those 
of the oscillators 1 and 2 because the Hamiltonian separates. Thus we have 

|0m,k) 2 = |0tn,l) <8> |0m,2) . (Bl) 

The one- mode squeezed states are governed by the same parameter z/2 : 



n=0 

The same expression holds for the ket \0in, 2). 

The overlap of this one-mode squeezed state with a one-mode coherent state is 

(v, l|0m, 1) = -= exp f-|fi| 2 /2 + zv\ 2 fi\ . (B3) 
y\at\ ^ ' 

According to [5], when taking small interactions into account, the density matrix of a 
one-mode squeezed state will preferable decohere into the mixture 

/j2 
-^Pi(t>i)h>i>(vi|. (B4) 

where the statistical weight is given by the probability to find a coherent state starting 
from the in vacuum : 

~Pi( v i) = \( v , l|0m, 1)| 2 = - — r exp (— \vi\ 2 + Re(zv{ 

\a\ V 



}_ p -2R\ r-ll^Ol^+Oin-V^Rih) 



Q: 



e --2Kl e -l{/2\a\ +U(n "'-Kill) _ 



In the second line we have introduced the real and imaginary parts of v\ in order to show 
that one gets an ellipse of great axis equal to \a\ which is oriented along the imaginary 
axis. The width of the small axis is 1/2, as in vacuum. 

For this decoherence procedure to be valid, as noticed in [12], it is important that the 
interactions do not break the coherence between k and — k modes, or equivalently do not 
mix and 02- 
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The product of two one-mode coherent states 1 and 2 is also the product of a coherent 
state for the k and — k modes: 



\vi)®\v 2 } = D 1 (v 1 )D 2 (v 2 )\0in,l) <g> |0m,2) , 

= D k (v)D_ k (w)\0in,k) <g> \0in, -k) , 
= \v, k) (g> \w, -k) , 

where the amplitudes are related by 



v = 



v\ + iv 2 



w = 



V\ — IV 2 



V2 



(B6) 



(B7) 



Finally, the product of the probabilities (B5) give the probability (32). Performing the 
change of variables from (v\, v 2 ) to (v, w) completes the proof. 

APPENDIX C: DIAGONALIZATION OF THE COVARIANCE MATRIX 



Introducing the position and momentum variables for each mode, i.e. = (q^ + 
ipk)/V^ and d_k = (q_k + iP-k)/V^, and defining the vector 



C f = (Ok Pk q-k P_ k ) 



one has the covariance matrix 



( n k + \ 


V Ci 







i 



C r 

^ n k + 2 
-c r n k + \ J 

2 



Ci \ 





(CI) 



(C2) 



where [, ] + is the anticommutator. Notice that (n k + \Y ~ \°k\ 2 > is a necessary 
condition for the matrix to have positive eigenvalues. 

The transformation Eq. (38) amounts to diagonahze this matrix: 

C = M l TM , 



(C3a) 
(C3b) 



The matrix T is the covariance matrix of the two thermal density matrices pth,a ® f>th,b 
in Eq. (38). The matrix M is the product of two local transformations and one global 
rotation R. The latter brings the covariance matrix C under a 2 x 2 bloc diagonal form. 
A product of local rotations Ri(0i) © R 2 {02) diagonahze each bloc, and the product of 
local squeezing Si(r{) © ^2(r2) brings the resulting matrix under the form T. Explicitly 
one has 



R(0) = 



( 


COS (p 





— sin (j) 





\ 







COS <j) 





— sin 


4> 




sin0 





COS 4> 







V 





sin^ 





COS (J 


* / 


( 


cos 9 
— sin( 


sin 6 
9 cos 6 


)■ 


S(r) 


= 



3 r/2 





-r/2 



(C4) 
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with the rotation angles d> = 7r/4, 6\ = 02 given by tan 29 = —Ci/c r , and the squeezing 
parameter r\ = —r2 defined by thr = — \ck\/(rik + 1/2). 

The eigenvalue of the thermal matrices is easily obtained by conservation of the 
determinant: 
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